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ABSTRACT 


Acoustic backscattering from a random rough water surface has 
been studied experimentally to test the application of two scatter- 
ing theories, a statistical description and a resonance approach. 
The rough surface was created by wind agitation over an anechoic 
tank. The wave height distribution was measured with a resistive 
probe and the wave slope distribution by optical glitter detection 
using a photocell. The distributions of backscattered sound pres- 
sures were recorded for surface roughnesses and sound frequencies 
corresponding to a very wide range of roughness conditions. Both 
statistical and resonance theories have regions of applicability and 
regions where they fail to predict the backscatter. Backscattering 
may be considered to be due to these two mechanisms, since a sum- 
mation of the predictions of the two theories fits the experimental 
data fairly well over the complete range of angles of incidence 
studied. A method of analysis is suggested to relate the doppler 


shifted continuous wave sound spectrum to the surface wave spectrum. 
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ACOUSTIC BACKSCATTER FROM Å 


RANDOM ROUGH WATER SURFACE. 
1. Introduction. 


Study of the limitations of acoustic propagation in the sea 
invariably reveals surface reverberation to be one of the chief 
noise sources. In recent years, many attempts have been made to 
produce a cogent theory to relate the scattered signal strength 
to parameters describing the sea surface, and these efforts have 
met with some success. However, results from measurements at 
sea show that no theory available at present can explain the 
scattering phenomenon completely. 

Two principal theories, a statistical theory and a resonance 
theory are examined briefly to show the predicted backscatter. The 
major part of this work is devoted to experimentally measuring 
the pulsed and continuous wave backscatter from a rough water 
surface and to measuring the necessary surface parameters as 
required by theory. The results are compared with the predictions 
of theory and with results of other experimenters both in the 
laboratory and at sea. 

Scattering theory has been developed largely with respect to 
electromagnetic radiation, yet the problem of acoustic scattering 
is completely analogous. Except for the phenomenon of polarization 
inherent in electromagnetic radiation, the reverberation problem 
in underwater acoustics is identical with the "sea-return" problem 
encountered in the operation of pulsed radar.at low angles over a 


rough water surface. 


1. 


Historical 


Several works of recent years have pointed out the nature of 
the backscatter phenomenon. In 1956, Urick and Hoover (1) conducted 
experiments at sea near Key West, Florida under various conditions 
of wind speed and sea state. Their data brought out essential 
qualitative details of acoustic backscatter. They noted that at 
grazing incidence, (transducer beam nearly horizontal) backscatter 
increased as the surface became rougher. At normal incidence this 
effect was reversed, the backscattered energy decreased as the rough- 
ness increased. In addition, they noted the existence of a precursor 
in their echo traces at normal incidence which led them to predict 
the existence of a subsurface volume scattering layer. Measurement 
of backscatter by Garrison, Murphy and Potter (2) in 1959 in Puget 
Sound, confirmed the general shape of curves of backscattering 
coefficient vs. grazing angle, although the latter used wind speed 
rather than surface roughness as a parameter. (That they are 
related parameters is beyond question.) It is noted that the latter 
paper used Urick's and Hoover's postulated volume scattering layer 
as a possibly explanation of poor definition of curves obtained at 
low wind speeds. In addition, the effect of the curvature of 
acoustic paths was considered as a significant source of error for 
the long range method used in the experiment. 

In 1961, Chapman and Harris (3), in deepwater north of Bermuda, 
extended the measurements at sea to the low frequency range using 
explosive sound sources. The backscattered echoes were recorded 
and played back through octave band filters. The curves of back- 


scattered signal strength vs. grazing angle for all frequencies 
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showed dependence on grazing angle at large grazing angles and 
virtual independence of grazing angle at small grazing angles. 
This effect they, too, attributed to the existence of a volume 
scattering layer. Since this region of independence of grazing 
angle was seen to be at a higher level for higher sea states, 
the possibility of the volume scatterers being bubbles was raised, 
although it was noted that, in this respect, the data were 
inconsistent. The existence of the volume scattering layer was 
deduced from the return of echoes prior to the expected reverba- 
tion from the shock. This pre-cursor had a return time which 
demonstrated a diurnal variation suggesting that the scattering 
layer might also be biological in origin. 

Scattering theory was developed in this period along two lines. 
Eckart (4) developed relationships between the scattering coefficient 
and the surface statistics of wave slope and wave height. Beckmann 
and Spizzichino (5) provided an excellent bibliography and summary 
as they examined the electromagnetic scattering phenomenon and 
related the scattered intensity, once again, to the surface para- 
meters of wave slope and wave height. In the work of Beckmann and 
Spizzichino and Shulkin and Shaffer (6), the significance of a 
surface roughness parameter such as the Rayleigh roughness criterion 
is noted, and the suggestion is made that correlation of backscattering 
with wind speed may be observed for fully developed seas but not 
necessarily for other sea conditions. 

In 1964, Clay and Medwin (7) examined the thesis that the volume 
scatterer might be bubbles and they calculated the bubble density 


required to produce such an effect. Barnhouse, Stoffel and Zimdar (8), 


15. 


in 1964 conducted experiments to show that bubble populations greater 
than that necessary to provide volume scattering did exist in 
Monterey Bay, California. 

Ån entirely different approach from the postulation of statistic- 
ally distributed scattering surface facets is taken by Wetzel (9). 
As with Marsh (10), he deals with the relation between sea scatter 
and the ocean wave spectrum. He raises the question of a resonant 
interaction of electromagnetic waves with the moving water waves 
to show Doppler shifts. The theory is derived only for the scatter- 
ing problem for relatively smooth surfaces. 

The principal strength of the statistical theory appears to be 
for large grazing angles and for specular scattering (scattering 
in the mirror direction). However, the theory fails to predict 
the backscatter observed at or near grazing incidence. The theory 
that the acoustical scattering at grazing incidence is not a surface 
scattering but a volume scattering phenomenon begs the question as 
to a postulated mechanism since the same phenomenon occurs in the 
scattering of electromagnetic radiation from the same surface. Hence 
it appears reasonable to demand that theory contain a prediction of 
an inflection in the form of the curve of backscattering coefficient 
vs. grazing angle to explain the scattering at grazing incidence. It 
must be borne in mind that the theoretical curves used for comparison 
are asymptotic forms for "smooth" surfaces and for "rough" surfaces. 
Even so, the departure of experimental evidence from the asymptotic 
form for rough surfaces raises doubts as to the validity of the 
statistical model. Fante (11) has shown how curves of the correct 


shape can be found by modifying the assumption concerning the form 
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of the correlation function for the surface. Beckmann and Spizzichino 
have assumed a Gaussian correlation function because it is consistent 
with the known Gaussian slope distribution which accompanies the 
Gaussian height distribution. Fante has pointed out that an expo- 
nential correlation function gives the correct shape for the curve. 
However, the actual levels of backscattering predicted by Fante's 
modification of Becktann's assumptions in a statistical scattering 
model appear to be seriously in error. The use of an exponential or 
Gaussian correlation function is discussed by Fung and Moore (19). 
The resonance theory fails at nearly normal incidence, although 
the backscattering is predicted, as measured, at grazing incidence 


to be nearly independent of the grazing angle. 
Summary of Experiment 


It appeared from the summary of scattering theory that there was 
room for a number of experiments under controlled conditions to show 
how the several theories and hypotheses fit (or fail to fit) exper- 
imental data. This study is a consideration of backscatter only 
although the theory makes predictions for acoustic scattering at any 
angle. The backscattering problem, then, is confined to the situation 
where the receiving hydrophone is identical with or contiguous to the 
projector. 

The statistical theory was tested by using pulse techniques. As 
seen from the theory, the backscattering coefficient is a function of 
the surface statistics and sound frequency through the definition of 
a roughness parameter. The surface statistics are the root mean 
Square wave height, the root mean square slope and the correlation 


distance. The three statistics are not independent and in this 
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experíment the wave height and wave slope only are measured. The 
wave height is measured by a resistance probe and the wave slope 

by optical methods similar to the Cox and Munk (12) and Schooley (13) 
glitter photography methods, Sound frequencies of 20, 50, 70, 100 
and 450 kHz were chosen so that the roughness factors approximated 
those experienced for a wide range of sea states for audio fre- 
quencies commonly used in underwater sound instrumentation. 

The resonance theory was tested using continuous wave tech- 
niques. Differential amplifiers were employed to cancel out the 
effect of the signal from the projector directly incident upon the 
receiving hydrophone. The resulting signal was processed through 


a fixed bandwidth filter to study the backscatter spectrum. 
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2. Backscattering Theory. 


At the present time the whole field of backscattering is in a 
State of uncertainty. Several theories have been postulated, but so 
far none has been sufficient to explain the observed phenomena for 
the full range of grazing angles. Two theories are considered in 
this thesis, the statistical approach as described by P. Beckmann 
and A. Spizzichino (5), and the "resonance" approach as used by 
H. W. Marsh (10), and L. Wetzel (9). A brief review of each is 
presented concluding with the expression for the Backscatter Factor. 
Å comparison of the theoretical and experimentally observed results 


is made in Section 4. 
Statistical Theory 


In this approach to the scattering theory the random rough 
surface is described by the statistical distribution of its wave 
heights about a mean level and by a correlation function that 
relates the wave heights at the various points on the surface. 
The distribution of wave heights is assumed to be Gaussian with 
mean zero and standard deviation O which is identically the root 
mean square wave height. The correlation function is here assumed 
to be 

C(t) = exp (-T2 72) 
where the parameter T is the correlation distance of the surface 
and T is the distance between any two points whose correlation is 


desired, 


19% 


The development of the statistical theory to be followed is a 
bríef summary of the derivation given by Beckmann and Spizzichino (5), 
Chapters 3 and 5, for electromagnetic radiation. Acoustic scattering 
theory is completely analogous to the electromagnetic, the pressure 
field p replacing the appropriate scalar component of the electric 
field E. 

It is convenient at this point to introduce the notation and 
geometry (Figs. 2.1 and 2.2) of the general scattering problem. 

The rough surface height is given by the function § = §(x,y). The 
mean level of the rough surface is the plane z = 0. Cartesian 
coordinates x,y,z, with origin O and unit vectors 1.,3 4 

used. The subscript 1 will denote all quantities associated with 
the incident field and the subscript 2 all those associated with the 
scattered field. 

6, = angle of incidence 


1 


93 = angle of reflection 


r = xíi + yl + 5(x,y)T, distance from origin of surface 
1 3 
coordinates to scattering element 


r) = distance from source to origin of surface coordinates 
ri - distance from source to scattering element 
r, = distance from orígin of surface coordinates to receiver 


r = distance from scattering element to receiver 
k, = kk, propagation vector օք incident ray 
k, = kk, propagation vector of reflected ray 
շը 
k = x propagation number 


P, = incident pressure field at scattering element 


pw scattered pressure field at hydrophone 
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Fig, 2.2 Scattering Geometry for Scattered Sound 
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We will assume Pj " Pio exp(ik, T) where 
P10 =፻ exp[i (kr; - Wt)] gives the amplitude and phase of the 
pressure field at the origin of the surface coordinates. The 
time dependence of the pressure field exp(-iwt) is suppressed. 
The pressure field P, at the hydrophone, assumed to be in the 


far field (i.e. the scattered rays are parallel), can be found 


by application of the Helmholtz integral. 


op 
za oe, =e 
Po ձո Ë Գո ۷ ծո | 95 (1) 
ջ 
Here p and — are the values of the pressure field and its 
5 on 3 


normal derivative at the surface. 





exp (ikr >) exp[i(kr, - k,r)] 
աջ» / T 
19 r, 


and fi is the inward normal from the surface. Making the substitution 


we arrive at 


3 c - ker)] exp[1(kr, - k.r)] 58, 
Po = Ar Ps an r " dn 


2 = 
Տ 
In applying the boundary conditions we make use of the Kirchhoff 


approximation that the plane surface conditions hold at each point 


and the assumption of a "pressure release" surface to get 


-1 exp (ikr, ) 
«այան 


= AN - e = 
ø - o S 


For simplification let v = ել - k, 


= A A 
v = k[ (s1n0,-sin0,cosQ)i, ^ (51ո0,51ոՓ)1շ - (6059) + cos0,)1,] 


22; 


I£ P10 is assumed constant over the ensonified surface and zero 


everywhere else we integrate and find 


-ip,gexP(ikr,) -— ub 
P» AL. x (k,n) exp(iver) dS (2) 
$ 


Equation (2) enables us to compute the sound pressure field P» 
once the random rough surface is determined since r and fi are 
functions of it. 


If the surface is smooth [£(x,y) = 0] 


աա ^ 
kjen -kcos0, 


Ver = wx + vy; 


and if the ensonified area is assumed to be rectangular, the scat- 


tered pressure field becomes 











1858166959) 
P, = nr, exp(ikr,) sinc ሃእ. sinc ا‎ 
sin VL 
where sinc v L = and A = 2L x 2W 
ጁ v L 
For specular scattering (9, = 6), v 7 % = 0) from a smooth surface 
the pressure field P20 at Էշ 18 
4 188 16 A ८050] 
Poo و‎ 2ፔ ՀՔ 07 (3) 
2 


Application of equatíon (2) to a random rough surface leads 


to the simple formula 


PA 3 
Po = Poo A exp(iver) dxdy (4) 


1 + cos0,cosO0, - sinO,sinO,cosqQ 
where F = վ - : 2 
5059) (6059) + 50893) 


25» 


This does not include the "edge effect" term which is 
negligible when the assumption A»»A^ is made. 

Since թշ is a complex quantity the mean square scattered 
pressure <P,P,*> must be calculated before we can talk about 
the scattered intensity. The brackets< > denote the mean value 
and the asterisk denotes the complex conjugate quantity. Հթշթշ» 
is related to the mean and variance at the scattered pressure 


field թշ by the relation 
<p,P, V[p,] +<p,>-<p, 


The mean pressure «p? scattered from a rough surface is now 


determined by taking the space average of equation (4). 


20] 7 
<p,> x —— exp[i(v.x + v y)] exp (iv 5) dxd 


Since we have assumed a Gaussian distribution for surface wave 


heights the mean value of exp(iv, 5) is 
A 
(exp (iv, £)/ = exp( - 2 ) 





= exp (- 5 ) 

whexe we define the roughness parameter 
g = of Ë (665591 + 60891. 

The mean scattered pressure then becomes 


> = - å W 
“P; p. F exp ( 2) sinc vL sinc K 


20 
This is the coherent component of the specularly scattered pressure 


field since it is in phase with the smooth reflected pressure "20 
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The variance of the acoustic pressure field թշ is now needed 
so it can be added to the squared mean pressure <p,>* <p> to 


complete the picture of the total scattered intensity. 
= Ի» - >. * 
Մ Լթշ) = <p,p,>> - <p,>- <p> 


Substituting P. and թշ" from equation (4) leads to 


2 f 
V 878 መ. P20 È) ||| E [iv Gu - x) բ iya = yl] (6) 


- >-< >. - >] 
[<exp(iv, (2, 25)) exp(iv, z,) <exp( iv z,)>|dx, dx dy, dy, 


This integration is performed by assuming the surface roughness is 
isotropic, changing to polar coordinates, and using the series 


expansion for the correlation function 


C(t) = ا‎ ITS) in evaluating Xexp[ iv, (z, - 22)». 


<exp[iv (z; = 22)» exp[-g(1 - C)] 


> 2298 
exp(-g) É exp(-m [T ) 


The result 


gro * ሮክ ° t | 
< mE" $ ` |Խշ 2 " ) & -at /T? 
ሃ[83] ہے ۔‎ 20 | I < ኘ Sam, E g) ) m exp| -mr /T ] | Tdr 
ծ m=] 
is integrated using equ. 14, Appendix B, Beckmann (5). The variance 
of the scattered pressure is then 
2 


TF T^exp(-g)P^q E ሯ 2 2 ×7 
VL», = A ES mm! exp[ - (v ሦ ተ+ Մ y) I / 4m] 
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Finally the mean squared scattered acoustic pressure is 


<pop}> » Հթշ»:Հթ» + ՄԼթշ1 


2 47 ራ 2 W 
*> up = 
MIL Poof exp(-g) | sine v,L sinc P + (7) 


2 = m 
ዝ፲ 2 2 på 
ia > 2 expl - (v, 1 ) 75] 


m=] 
The first term is the coherent component of the scattering and is 
important only in the specular direction and, there only if g < 1. 
The second term, the incoherent component, is significant for all 
nonspecular scatter including backscatter. 


For scattering in the specular direction, (6, = 93 and ọ = 0), 


<p,p2> 
the ratio neo is important. Great simplification of 
P 
20 


equation (7) occurs since F = 1, and sinc Մե = sinc E = l. This 


leads to 





<p, p$ > ብር » m 
2 " exp(-g) +- exp(-g) سس‎ 


The bracketed portion of the second term in the above equation 
is defined as S(g). This factor takes on the extreme values of 


S(g) =e for g<.1 


۷ 


S(g) = 1/g for ४ 210 


for the relatively smooth and very rough surface respectively. 


For the smooth surface (g<.1) 


27 
meg 
اکا‎ 
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For the rough surface (g>10) 


* 
“8283 ብ 3 ուշ 


E Ag 





| P20 


or substituting for g and putting in terms of the mean square slope 


we have 


*> 
ጨበር ካህ X 


| IP, ol f 812 260520 


For the backscatter case v? + Y = 4k sin^e, and 


F = (८059)) since 6, = -90, and oç = 0. The surface backscattering 


Ji 2 
factor is defined 


ՀԹշթշ- 


T ————— छी 


2 
LA 


H 
>| N 


5, (8) 


We can now find the equations for the backscattering factor Շր) 
for a relatively smooth surface (g<<l) and a rough surface (g>>1). 
For the relatively smooth case (g««1) the first (coherent) 


term in the brackets of equation (7) is negligible for 0, # 0; 


(it must be considered when 91 = 0). 


the first term of the series need be taken hence 


Furthermore since g<1 only 





på F^nT^g 

20 2 2 nå 

<P, P$? "ium EE exp[ - (v ጅ n) /^ - g] 
lón^r^o 299 0 2 2 

<P,P$? काठ 2 exp(-k T2sin 9; - ջ) 
እ ፍዐ5 9 À 


ZUR 


The backscattering factor S. then becomes 


B 
35222 
TG 
Տ A TONENE ھی‎ kai sin O - g) (9) 
B 4 1 


For the rough surface (g>>1) direct application of 
equation (7) is not reasonable because the second term 
converges too slowly to be of practical use. We must return 


to equation (6) and make the approximation 
4 2 -- 
<exp[ iv, (z -z,) ]> = exp[-g(1-C)]J  expl-g /T ] 


where C= 1- T/T. 

Using the same method of integration as before and noticing that 
the second term in the brackets (equivalent to exp[-g]) approaches 
zero for g>>l we neglect it and arrive at 


jas 
TF 0 


of: 
k d بے رید‎ 
ብ Ea expl-(v, * P2 /4g] 


For backscattering this is equivalent to <p,på> since the coherent 
term is negligible. Substituting in equ. (8) the backscattering 


factor becomes 
Տը = (BTA ^cos 9.) | exp[ -tan^0, /2 uem (10) 


where ነዎ is the mean square slope. 
In general, the data taken from backscatter experiments at sea 
appear to follow equation 10, decreasing rapidly with increasing 9 
at angles of incidence up to 40° measured from the normal. The 
equation does not however, predict the virtual independence of 8 


1 


that has been observed for angles of incidence greater than this. 
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Resonance Theory 


In the resonance approach to backscatter the random rough 
surface is described by a wave energy spectrum. This spectral 
function E(f) is proportional to the surface wave "energy" per 
cycle bandwidth (em /aec). The nature of the sea surface 
spectrum (Fig. 2.3) shows that there is very little energy in 
the band of frequencies below the characteristic frequency f ° 
determined by wind speeds; that most of the wave energy lies in 
a small band of frequencies centered around f and that at higher 
frequencies the curves fall as f, where a common value for n has 
been found to be 5, Phillips (18). 

Wetzel (9), in his approach to the problem, postulates that for 
the relatively smooth (g<<l) surface the major backscatter radiation 
does not come from the incoherent (scattered radiation arriving at 
receiver with uniformly distributed phase) component caused by the 
randomness of the scattering surface. Rather he suggests the 
backscattering is due to a resonant interaction between the incident 
radiation and some component in the water-wave surface. This inter- 
action leads to a constructive interference which is physically 
similar to Bragg diffraction, and which, because of the surface wave 
velocity, manifests itself in a coherent doppler shifted frequency 
in the non-specular directions. 

Using Bragg's Law 

2LsinO = nÀ n 122599 0 e 
one can compute the appropriate surface wave length for a given 
frequency of acoustic radiation and a known angle of incidence 9]: 


For the first order effect (n = 1) the phase velocity for this 
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Fig. 2.4 Resonance Interaction 


E 
ዘ 


particular wavelength is then computed (taking into account surface 
tension if the wavelength is less than a few centimeters, Fig. 4.29) 
and the doppler shift calculated from the relation 


Af = 2Vsin@, f/c (11) 


where V is the phase velocity of the surface wave. 

Wetzel (9) derives an expression for the backscatter cross 
section per unit area, using the Kirchhoff approximations, and based 
on Chapter 3 of Beckmann and Spizzichino (5) modified slightly to 
include a term that Beckmann had identified as negligible. With this 
modification he arrives at the following expression for the back- 


scatter cross section per unit area (Wetzel (9), equation 41) 
d "ES 
ճաշ 


o, /À - k'g K. ^ E [£()] (12) 


where K. = 2ksin9, is the resonance wave number, À is the surface 


1 
area ensonified, g is the acceleration of gravity, and E Լ 08.2) 


is the value of the surface wave spectral energy at the resonance 


wave number. Substituting the values for K. and g leads to 
2 3 
هير‎ = 11.07k : E [£(K ) )/sin 281 cgs units 
Using the definition of the backscattering factor (8) it can be 


shown that this is equivalent to o, ,/^ = ձոՏը or 


572 
.881k 


i sin? ^g. 


The level of backscattered sound predicted by equation (13) 


. E [፻(፳ ) ] cgs units (13) 


is examined quantitatively in Section 4 by comparing with the 
observed values. The values for E L£(K )] are obtained directly 
from the wave energy spectrum measured (Fig. 3.33) for the surface 


by Cohen and Scheible (17). 
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In general Wetzel's equation (41) is a plausible explanation 
of the low level backscatter that has been widely observed for 
small grazing angles. However, at normal and near normal incidence 
it fails to give a reasonable prediction. The theoretical expression 
for the backscattering factor from the resonance theory appears to 
have a singularity at the normal incidence prediction but this is 
removed by the nature of the energy density in the surface wave 
Spectrum which has no energy in the spectrum for a surface wave 


with zero frequency. 
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3. The Model Sea Surface 


Generation of Rough Surface 


All experiments were conducted in a fresh water filled anechoic 
tank 25' x 6' x 7' deep. The rough surface was produced by a combina- 
tion of up to three centrifugal blowers. The tank surface was covered 
for a distance of ten feet from the point where the fan discharge 
struck the water surface vertically, by polyethylene sheets held in 
place by a slotted angle framework and plywood to form a plenum 
with its top 3% inches off the surface of the water. This provided 
a short fetch for the wind to develop a rough surface. When it 
seemed necessary, a wire mesh "beach" was inserted one foot from the 
end wall of the tank to suppress standing waves in the water surface 
due to the small reflection from the end of the tank. The arrange- 
ments are shown in Figs. 3.1 and 3.2. Three Locations or stations 
are marked in Fig. 3.1 to indicate the centers of areas of scattering 


used in the conduct of the experiments. 


Wave Height Statistics 


The statistical theory assumed that the rough water surface 
is a stationary random process and that the distribution of wave 
heights is Gaussian with zero mean and standard deviation, O. It 
is desired to test this assumption from a theoretical point of view 
and then to test its validity experimentally by measuring the actual 
wave height distribution of several rough water surfaces. 
Theory 

For a simple sinusoidal surface in one direction, it can be 


shown that the probability density function for the height is 
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bimodal and of the form 


P(£) = ےس چمےہ‎ for E<A 
à TT ል" - हु“ » 


where Å is the amplitude of the sinusoidal surface. The distribution 
is shown in Fig. 3.3, and is anything but Gaussian in nature. How- 
ever, the evidence of our eyes and the study of wave spectra show 

that the surface of the sea is not å simple one dimensional sinusoidal 
disturbance. If we consider the height at any point to be due to the 
superposition of Յ number of wave motions, independent in space and 
time, the application of the central limit theorem shows that the 
values of wave height tend toward a Gaussian distribution. Cox and 


Munk (12) show the development of just such a distribution. 
Experimental 


Since the theory depends not on the mean wave height but only on 
the root mean square wave height, the stationary surface level is 
taken as zero and the variance of the height distribution is measured. 
Figs. 3.4 and 3.5 show the apparatus used to measure the distribution. 

Probe. The probe is a piece of tinned copper wire, 0.80 mm in 
diameter and 25 cm in length, held with approximately 17 cm of its 
length below the surface of the water. The small diameter is necessary 
to minimize the formation of a wake behind the probe due to motion of 
the surface layer of water. The return is provided by a second 
electrode 3.14 mm in diameter and 45 cm in length, which is completely 
submerged at all times. This second electrode is connected to ground 
and in fact is parallelled by many other ground connections in the 
plumbing of the tank, so that its position and size do not affect 
the operation of the device to any extent. The impedance of the 


cell formed by the short surface probe and the ground return is 
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مابس 


- — ah 


then a function օք the amount of the probe below the surface and 
hence is a function of the wave height at the probe. The physics 
of the operation of the probe have not been thoroughly examined, 
although it is clear that conductive effects rather than capacitive 
effects predominate in its operation. Farmer and Ketchum (14) have 
given 9 theoretical derivation of the sensitivity of a device using 
similar methods, though of a much larger scale. What was clearly 
necessary, however, was an accurate calibration to determine the 
linearity of the device in measuring wave height and a frequency 
response test to determine whether the varying height of the water 
surface could be accurately changed into an electrical analog. 

Electronics. Fig. 3.6 shows a block diagram of the system 
used in the determination of the wave height distribution. The 
device makes use of the output impedance of the oscillator by 
causing the output to be shunted by the variable impedance of the 
probe. This varies the output voltage of the oscillator. The 
resultant signal is an amplitude modulated 10 Hz carrier where the 
modulation contains the wave height information. The AC carrier is 
necessary to prevent the electrolytic polarization of the probe which 
would occur if DC were used, 

The detector is a simple full wave rectifier followed by a T 
section RC filter as shown in Fig. 3.7. The filter must have a low 
pass characteristic in order to remove the carrier frequencies but 
not such a low cut-off as to cause undue phase distortion of the 
modulation signal. The differential amplifier is required to remove 
the large DC portion of the signal and to increase the signal level 


for the pulse height analyzer. 
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Oscillator 

Probe 

Detector 

Differential Amplifier 
Variable DC Source 
Sampler 

Pulser 

Pre-Amplifier 

Pulse Height Analyzer 
Tape Recorder 


Wave Analyzer 


HP 650A Test Oscillator 

see Fig. 3.5 

see Fig. 3.7 

Tektronix 1A7 in Tek 545 Oscilloscope 

HP 467A Power Amplifier 

see Fig. 3:8 

GR 1217B Unit Pulse Generator 

HP 467A Power Amplifier 

RCLiac 128 Scaler Analyzer 

PI-6200 Portable Instrumentation Tape Recorder 


GR 1900A Wave Analyzer 


Fig. 3.6(a) Equipment List for 
Wave Height Measurement System 
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Fige 3-8 Sampler Unit 
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The sampler consists of a field effect transistor which is 
controlled by a high speed silicon switching transistor, as shown 
in Fig. 3.8. The pinch-off voltage of the FET is removed when the 
pulse causes the switching transistor to saturate thus providing 
a very short rise time at the opening of the sampling gate. This 
fast rise time is required for the input to the pulse height 
analyzer. 

The pulse height analyzer is a RCLiac 128 Channel Scaler Analyzer 
designed principally for nuclear physics applications. The scintilla- 
tion crystal, photomultiplier and preamplifier section were removed 
and the input provided to the variable gain section through a separate 
isolation preamplifier. The analyzer is digital, providing a cathode 
ray tube output of the pulse height distribution both pictorially, as 
a frequency diagram or histogram, and as a decimal numeric readout 
for each channel. The analyzer detects the rise of a pulse and 1% 
microseconds later it picks off and holds the input voltage. Simul- 
taneously an internal 2 MHz counter is started and a rising ramp 
voltage is generated. When the ramp voltage exceeds the held input 
voltage, the count is stopped and a single entry is made into the 
cell in the ferrite core memory whose address is the count at stop 
time. The input is then re-enabled for the next pulse. The system 
dead time is 70 + n/2 microseconds, where n is the address of the 
memory core for the count. It is seen that the input pulses must 
have a sufficiently small rise time to rise to the full value in 
less than 1% microseconds. In addition, the analyzer must be fed 
by some pulsing scheme other than allowing it to sample whenever 
the input is enabled, otherwise the sampling rate would not be 


independent of the pulse heights and the lower pulses would be 
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Fig. 3.9 Surface Probe Characteristic 





favored. In this case, the sampling frequency for the wave heights 
was set at 1 kHz which is well beyond the limit dictated by the 
Shannon sampling theorem to prevent aliasing in detecting the 
highest frequencies in the wave height signal. Since the highest 
n is 128, the analyzer is limited to a maximum regular sampling 


frequency of 7.4 kHz. 


Gali ation. Examination of the characteristics of the output 
voltage of the wave height measuring device under slowly varying 
conditions showed that the output voltage vs water height was of the 
form shown in Fig. 3.9. It is clear that the relationship is any- 
thing but linear but over a small range such as the two to three 
centimeter variation expected in the model rough surface, the 
linearity is very nearly preserved. It is also noted that linearity 
improves as the immersed length of probe increases, but only at the 
sacrifice of voltage sensitivity of the device. Hence a compromise 
was made at an immersed depth of 17 cm at which the sensitivity of 
the probe with the supplied voltage at 2.0 volts was 1.3 mv/cm. 

It is not sufficient to use a static calibration of the device, 
however, unless it is clear that the device will operate accurately 
over the range of frequencies in the surface wave spectrum. Since 
it seemed more practical to move the probe up and down in the water 
in a sinusoidal motion rather than to cause the water to move past 
the probe in a sinusoidal motion, the former course was chosen. 

Fig. 3.10 shows a diagramatic sketch of the apparatus used. The 
probe was driven vertically by a simple cam mechanism operated by a 
variable speed motor. The amplitude was fixed at about 3 cm by the 


cam arrangement but the frequency was continuously variable up to 
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Fig. 2.10 Dynamic Calibration of Wave Height heasuring Device 
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8 Hz. The actual position of the probe was determined by a linear 
potentiometer supplied by a constant voltage DC supply. The 
potentiometer output was led to the X-input and the detected wave 
height led to the Y-input of a cathode ray oscilloscope. The 
resultant 1:1 Lissajou pattern was analyzed at several frequencies 
for linearity, sensitivity and phase distortion. It was expected 
from the static calibration that the linearity would be good, but 
the effects of surface tension to cause the inevitable meniscus about 
the probe causing a lagging response, were unknowns. 

The results of the first calibration showed an unacceptably 
large phase lag at any frequency above 3 Hz. Since this phase lag 
could be cauSed by the characteristics of the filter in the detector 
as well as by the effects of surface tension it was decided to attempt 
to modify the detector. Fortunately it was found that simply reducing 
the time constants of the filter to the values determined by the 
parameters in Fig. 3.7, forced the poles of the transfer function of 
the system far enough out along the real line in the complex frequency 
plane to provide an acceptably low level of phase distortion. It 
was evident that the poles of the filter system dominated the pole 
introduced by the surface tension. After adjusting the filter 
parameters it was determined that at 8 Hz, the highest frequency at 
which calibration with the mechanical oscillator was possible due to 
limitations of the controller, the phase lag was 3 degrees. Since 
the system is essentially third order, the response can be extrap- 
olated on a worst case basis on the assumption that all three roots 
of the system equation are identical. Fig. 3.11 shows such a 
construction as a Bode diagram to show the worst response in 


magnitude and in phase. It is seen that the response of the device 
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is accurate to less than 1 db at frequencies up to 100 Hz at the 
very worst and is probably much better. It should be noted that 
the filter characteristics decided upon to optimize the frequency 
response of the device were not the best for suppression of the 
20 kHz ripple in the detected signal. However, the change in the 
standard deviation of the wave height distribution due to the 
ripple is quite small (1.4 mv) and can be corrected using the 
summation principle of variances for the sum of random variables 
which are uncorrelated. 

The conclusion of this dynamic calibration procedure is that 
the response at frequencies expected in the surface wave spectrum 
is flat, without phase distortion, and that it is sufficient to 
calibrate the device under static conditions and to apply that 
sensitivity to the entire spectrum. The pulse repetition 
frequency of the sampler was 1 kHz and the sampler pulse length 
set at 25 microseconds. The device was calibrated statically in 
smooth water by adjusting the immersed depth of the probe by turning 
the 20 pitch adjusting screw. The attached scale, which can be 
seen in Fig. 3.5, indicated the number of turns of the screw and 
the head of the screw was indexed to ensure the accurate setting of 
an integral number of turns. A typical calibration curve is shown 
in Fig. 3.12. It will be noted that the calibration yields a 
sensitivity in channels/cm, an odd unit perhaps, but convenient 
in the processing of wave height distribution data from the pulse 


height analyzer. 


। > 


49. 


१०6 


Wave Height 


-.6 


50 100 
Channel Identification 


Fig. 2.12 Wave Height Measurement Calibration 
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Results and Conclusions 


The wave height distribution was studied at three stations 
in the tank (see Fig. 3.1). The process is assumed random, 
homogeneous and stationary, that is the height distribution statis- 
tics are constant in time and over the area ensonified in the 
scattering experiment. It is further assumed that the statistics 
of the time distribution are the same as for the spatial distribution. 
Measurements of the samples of wave heights were taken for 197 seconds 
giving a distribution of 2.0 x 10° samples of wave height in that 
period. In each of the stations, the wave height was measured with 
combinations of one, two and three blowers operating. 

Capillaries were evident at station 1 with two and three blowers 
operating, but these had died out almost completely by the time the 
wave motion had reached station 2, eight feet further downwind. 

Figs. 3.13 to 3.18 show the height distributions at each of the 

three stations. The measured distributions are compared with Gaussian 
curves with the same (zero) mean and standard deviation, normalized 
to have unit area. Table 3.1 summarizes the measured wave height at 
each scattering location (stations 1 and 2) as well as at station 3 
which was the site of the continuous wave scattering experiment. The 
RMS wave height was calculated from the distribution data by the 
general purpose computer (CDC 1604 - FORTRAN) using the sensitivity 
calculated from the static calibration curve taken at the time and 
station where the wave height was recorded. A typical polaroid 


photograph of the pulse height analyzer readout is shown in Fig. 3.19. 
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TABLE 3.1 


RMS WAVE HEIGHT - CM 


1 Blower 2 Blowers 3 Blowers 
Station 1 0.099 0.121 0.206 
Station 2 - 0.129 0.216 
Station 3 - - 0.362 


The assumption in theory that the wave height distribution is 
Gaussian is very nearly correct for the wind agitated water surface 


in the tank. 
Wave Slope Statistics 
Theory 


In describing a random rough sea surface it is often more conven- 
ient to work with the slopes of the surface rather than with the 
heights above a mean level. In Appendix D of the book by Beckmann 
and Spizzichino ©) it is shown that a surface whose waveheights are 
normally distributed has its wave slopes normally distributed if the 
correlation function is Gaussian. More specifically, if the correla- 
tion function is C(T) = >= ም) and the wave height distribution 
has zero mean and standard deviation O the mean square slope turns 


out to be 


where नि, is the RMS slope angle. 
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Fig. 3.20 Apparatus for Measurment of Wave Slope Statistics 
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Experimental 


The method used in obtaining the statistical distribution of 
the water surface slopes was a modification of a simple optical 
method described by A. H. Schooley (13). Rather than obtain the 
Statistical distribution of the wave slopes from a series of 
"glitter" photographs (an extremely time consuming method), the 
following technique, as illustrated in Fig. 3.20 , was adopted. 

À simple box camera was constructed *d a Hoffman type 55C 
(.2" x .2") photocell placed at the center of the film plate position. 
A 100 watt zirconium arc point source of light was placed next to the 
lens of the camera (Figs. 3.21 and 3.22) and was assumed to be 
coincident with it. A convex lens was placed approximately 12" in 
front of the point source in order to concentrate the light intensity 
in a spot of diameter approximately 20 centimeters on the water 
surface. Provision was made so that the camera axis, passing 
through the lens center and the photocell, could be oriented to any 
angle from 0 to 20 degrees with the vertical. The whole apparatus 
was mounted on a movable platform over the tank with the camera lens 
at a height of 62" above the water. Measurement of the water surface 
slope was accomplished by adjusting the elevation rod so that the 
camera axis was at the slope angle desired. The inclined camera 
was then positioned horizontally so the measurement could be taken at 
the center of the scattering area (defined in Fig. 3.1). 

Every time a wave facet of the roughened surface was perpendicular 
to the axis of the camera (i.e. set at the slope angle desired) light 
from the point source was reflected to the photocell. The voltage of 


the cell, being proportional to both the intensity of light falling 
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on it and the partial area illuminated, varied with the movement of 
the water surface. The signal was essentially zero when the cell 
was dark, rising and falling rapidly (Fig. 3.25) every time a 
sparkle of light was received from the surface of the water. The 
signal was recorded directly, and after integration, for a two 
minute period to obtain the average value. This value was then 
plotted on a histogram for wave slope angles. The position of the 
camera was changed to a new slope angle and the process repeated 
until the complete histogram was obtained. 

The instantaneous signal from the photocell was approximately 
50 to 500 microvolts varying with the slope angle. This was 
amplified with a Sanborn D. Շ, ET i amplifier set with a 
gain of 1200 and then integrated with a G. A. Philbrick Operational 
Amplifier configured for integration (Figs. 3.23 and 3.24). The 
differential amplifier made this whole technique possible by 
reducing the noise level sufficiently to give a suitable signal- 
to-noise ratio. Both the instantaneous signal and its integrated 
value was recorded with a Brush Recorder oscillograph (Fig. 3.26). 

There was a bit more involved in the technique used to integrate 
the signal due to the small drift in the zero level of the output of 
the differential amplifier. This zero level drift was a major problem 
and several D. C. amplifiers were tried before one with a tolerable 
drift was found. After the differential amplifier was adjusted to 
minimize the D. C. component a small and variable D. C. component 
still existed. The signal from the photocell was superimposed on 
thís "noise" so it had to be corrected in the integration process. 
This was done by integrating the background level (point source off) 


for 30 seconds before and after each one minute integration of signal. 
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68. 


The signal plus background was integrated for two one minute periods 
and was then corrected for the background level integration. The 
average of these two corrected readings was plotted on the histogram 
for the slope angle. To find the accuracy of each point plotted on 
the histogram the standard deviation for several readings at a given 
angle was measured. These errors were found to vary from 4% to 9% 


and are noted on the histograms presented. 


Calibration 


The output voltage of the cell was assumed to be directly propor- 
tional to the portion that was illuminated (at constant light intensity). 
The plot of voltage output as a function of exposed area was found to 
be linear. 

The change of light intensity at the photocell due to increased 
spherical divergence loss for the increased range of the source and 
lens at large slope angles was measured. It was found that the 5% 
increase in range for the maximum slope angle resulted in a 2% 
decrease in voltage. This allowed the distance from the lens to the 
rough sea surface, and therefore the intensity, to be treated as a 
constant so no corrections had to be made to the photocell's output 
voltage. 

To obtain better slope angle resolution the photocell was 
partially masked with opaque masking tape. With the masked photo- 


cell the resolution for slope angle was + .2 degree. 


69. 


Results and Conclusions 


The slope angle distribution was taken in both the upwind- 
downwind direction and the crosswind direction for both scattering 
areas with two and three blowers on. The results are tabulated 


in Table 3.2. 


TABLE 3.2 MEAN SQUARE SLOPE 


Station I > > ep SG S | > 
2 blowers . 007 69 . 0028 .01007 3.25 
3 blowers . 00975 . 00680 . 01655 1.44 

Station II 


3 blowers . 00935 . 00364 ን ያየን 2.57 


The general results (Figs. 3.27 - 3.32) found for our model 
compare favorably with those measured at sea by Cox and Munk (12). 
The skewness in the downwind direction of the slope distribution 
reported by them was observed for our rough surface (Fig. 3.29). 

The crosswind distributions as well as the upwind-downwind distribution 
for the relatively smooth surface were all found to be approximately 
Gaussian as at sea. The ratios of the up-downwind to the crosswind 
component of mean square slope varied from 1.44 to 3.23 whereas those 
found by Cox and Munk (12) at sea ranged from 1.0 to 1.9. Our ratios 
were somewhat larger than at sea because of the channeling effect of 
our model. The mean square slope in our model corresponded to a wind 


speed of approximately 5 knots. 
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Fig. 2922 Energy Spectrum of Rough Surface 
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Surface Wave Spectrum 


The wave height signal from the wave height measurement system 
can be analyzed to determine the wave spectrum. This.description 
of the rough model sea surface allows comparison with ocean wave 
spectra for the purpose of assessing a scaling factor for the model 
surface and is also included for identification of the backscattered 
energy spectra measured in Section 4. 

This experiment was carried out as an advanced laboratory project 
by Cohen and Scheible (17) who used the identical system described 
in Fig. 3.6 except for the method of analysis. The detected modula- 
tion envelope was recorded at 0.375 ipa and then played back at 37.5 
ips into the wave analyzer with the fixed bandwidth selected at 3 Hz. 
In this way the spectrum was shifted to audio frequencies and expanded 
so that greater resolution within the spectrum is achieved. The 
resultant spectrum is shown in Fig. 3.33. It is noted that the energy 
is concentrated in a band of frequencies 6 Hz wide with a pronounced 
peak at 3.8 Hz and falling off at frequencies above the peak with an 


approximately ፲፪. dependence to 6 Hz and 5:7 beyond 6 Hz. 
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4, The Backscattered Radiation 


The statistical scattering theory predicts a relation between 
the backscattered intensity and the surface statistics and the 
frequency of the acoustic energy scattered. It also predicts the 
existence at normal incidence of a coherent and an incoherent 
component of backscattered sound. (At normal incidence, back- 
scatter and specular scatter are identical). For a relatively 
smooth (g<<l) surface, the distribution of the heights of the 
pulses scattered in the specular direction is predicted to be 
Gaussian. The distribution of energy scattered in other than the 
specular direction is predicted to be Rayleigh since the scattered 
energy is incoherent. For rough (g>>1) surfaces, the specular 
scattering is predominately incoherent and all distributions are 
expected to be nearly Rayleigh. The measurement of the scattered 
echoes has been designed to display the distribution of echo 
strengths and to determine the average backscattered itensity. 

The latter is calculated from the square of the mean plus the 


variance of the scattered pulse amplitude distribution. 
Experimental 


In the conduct of the experiment, it is extremely important 
to identify the received signal as backscattered energy from the 
particular angle considered. It is possible to measure noise, 
echoes from other portions of the tank or transducer support or 
from the side lobes or even off-axis in the main beam. Particular 
care must be taken with the geometry of the apparatus and with 


choice of electronics equipment to discriminate against all signals 


TOS 


which are not truly surface backscatter. With respect to electrical 
noise, low noise figure electronic apparatus is essential as is the 
use of bandpass filters to reduce the effects of wideband noise. 

So far as the geometry is concerned, any pulsed technique can 
discriminate against echoes from the tank walls or floor by the 

simple expedient of making the distance to the scattering area 

shorter than any other distance in the tank, or at least making the 
distance to the scattering area different from the distance to other 
possible targets so that the surface scattering and the spurious 
echoes do not coincide in time. However, in a model study of back- 
scatter, there is å spurious echo from the surface for which this 

is not always possible. At angles near normal incidence, the normal 
incidence reflection from off-axis in the main beam may be very close 
in time to the scattered echo, and the amplitude of this reflection 
may be larger than the scattered echo. To reduce the errors, the 
projector must have a very narrow beam pattern, yet, in a tank 
experiment this dictates that the dimensions of the scattering area 
will not be large with respect to the surface correlation distance. 
For a beam incident. on a relatively smooth surface the strength of 
the backscattered ا‎ near normal incidence is very sensitive to 
small changes in the angle of incidence. Since all other backscattered 
intensities at other angles of incidence are referenced to the 

reading at zero incidence, errors in the measurement of the scattering 
intensity at this sensitive spot will affect the level of all results. 
It is considered that this is 8 principal source of experimental error 
in this work. The level of electrical noise and spurious echoes was 
measured and these effects are indicated by a noise line on the 


curves of backscattering factor vs. angle of incidence (Figs. 4.17 to 
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Timing Oscillator 

Tone Burst Generator 

Signal Oscillator 
Plug in Unit 

Shunt Gate 

Power Amplifier 

SCR Gate 

Projector 

Delay Generator (2) 

Delayed Pulse Generator (2) 

Hydrophone 

Amplifier 

SKL Filter 

Receiver Blanker 

Detector 

Sampler 

Isolation Amplifier 


Pulse Height Analyzer 


HP 650A Test Oscillator 

GR 1396A Tone Burst Generator 

HP 3300A Function Generator 

HP 3302A Trigger/Phase Lock Module 
see Fig. 4.4 

GR 1233A Power Amplifier 

see Fig. 4.5 

EDO Model 327 Transducer 

GR 1217B Unit Pulser 

GR 1217B Unit Pulser 

Atlantic Research LC-32 

HP 463A Precision Amplifier 
Spencer Kennedy Model 302 Electronic 
see Fig. 4.6 

see Fig. 4.7 

see Fig. 3.8 

HP 467A Power Amplifier 


RCLíac 128 Scaler Analyzer 


Fig. 4.2(a) Equipment List for 


Transmitter and Receiver System - 70 kHz 
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Timing Oscillator 
Signal Oscillator 
Tone Burst Generator 
Shunt Gate 

Power Amplifier 

SCR Gate 

Projector 

Delay Generator (2) 
Delayed Pulse Generator (2) 
Hydrophone 
Pre-Amplifier 
Amplifier (2) 
Passive Filter 

SKL Filter 

Receiver Blanker 
Detector 

Sampler 

Isolation Amplifier 


Pulse Height Analyzer 


HP 650A Test Oscillator 

HP 650A Test Oscillator 

GR 1396A Tone Burst Generator 
see Fig. 4.4 

GR 1233A Power Amplifier 

see Fig. 4.5 


4 x 4 array BaTiO, special NPGS 


4 
GR 1217B Unit Pulser 
GR 1217B Unit Pulser 


3 x 3 array BaTiO, special NPGS 


4 
HP 463A Precision Amplifier 

HP 467A Power Amplifier 

Allison Model LC Variable 

Spencer Kennedy Model 302 Electronic 
see Fig. 4.6 

see Fig. 4.7 

see Fig. 3.8 


HP 467A Power Amplifier 


RCLiac 128 Scaler Analyzer 


Fig. 4.3(a) Equipment List for 


Transmitter and Receiver System - 450 kHz 


4.21). Unless these particular precautions are observed in any 
scattering experiment, the small grazing angle effects, which 
depart from theory, and which have been observed by several 


experimenters (1,2,3) may well be illusory. 


Electronics. Fig. 4.1 shows the electronic equipment used in 
this part of the experimental work. Figs. 4.2 and 4.3 are block 
diagrams of the systems used in measuring the backscattered radia- 
tion from the rough surface at 70 kHz and 450 kHz. They differ only 
in the use of additional amplifiers at the higher frequency to 
compensate for attenuation inherent in the sampling circuits at the 
higher frequency, and in the use of a different signal generator in 
each of the experiments. 

The tone burst generator used in the experiment at 70 kHz was a 
HP 3100A Function Generator with a 3302A Trigger/Phase Lock plug-in- 
unit. This provided 67 db separation between the tone burst level 
and the output noise. For the experiment at 450 kHz, a separate 
signal oscillator was used and the tone burst formed by a GR 1396A 
Tone Burst Generator. The GR tone burst generator provides only 44 
db separation between tone burst level and leak-through signal. 
Particularly in the case where the transducer is used both as a 
projector and a hydrophone, there must be at least 80 db separation 
between the tone burst level and the noise or leak-through signal. 
An additional 27 db was provided by a supplementary shunt gate, 
and in both cases, noise from the power amplifier during the off 
period was suppressed by use of a switch incorporating a silicon 


control rectifier. 


86. 


The supplementary shunt gate is shown in Fig. 4.4. It is 
a simple emitter follower followed by a shunt switch with provision 
for pedestal adjustment. The switching circuits of the GR Tone Burst 
Generator are very similar. With the transistors used, one stage 
provides 27 db separation between the signal levels measured in the 
on and off modes. The key factors in determining this ratio are 
the transistor resistances in the saturated and cut-off modes. 

The power amplifier used in the pulsed experiments was a GR 
1233A Power Amplifier. This unit provides a high gain and a wide 
frequency response with signal levels up to 120 volts peak-to-peak. 
This amplifier is subject to oscillation unless particular care 
is taken in location and layout of input and output circuitry to 
avoid feedback. It is also particularly sensitive to poor pedestal 
adjustment in the input signal. Any amount of pedestal voltage 
causes distortion of the output pulse. 

The Silicon Control Rectifier (SCR) gate is used at the higher 
power level to isolate the power amplifier when the system is in the 
listen mode. The low output impedance of the power amplifier is an 
advantage when supplying power to the transducer, but it is a decided 
liability if the transducer is being used as a hydrophone as well, 
since the signal is shunted by the power amplifier output impedance. 
Secondly, the noise output of the power amplifier with zero input 
signal is still sufficient to mask the scattered radiation and must 
be isolated. Fig. 4.5 shows the SCR gate. The gate is self operating 
and requires no external gating pulse. The diode bridge ensures that 
the SCR faces signals of the correct polarity for both positive and 
negative halfcycles of the tone burst. Initially the SCR is not 


conducting. The positive rise of the input signal biases the cathode 
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Fig. 4.6 Receiver Blanker Circuit 
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Fige 4.7 Detector Circuit 
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Fig. 4.8 Regulated Power Supply 
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gate and the SCR turns on. As the current falls after the first 
half-cycle, the critical hold-on current is passed and the device 
turns off. If the tone burst continues, the process is repeated 

for each half-cycle. If the tone burst is finished, the device 

stays turned off. The resistance in the on mode is very low 

(10 ohms) and in the off mode is very high (2 megohms). It 

should be noted that the gate voltage required to turn the SCR on 

is about 2 volts. For a large signal, this causes negligible 
distortion. For a small signal, the distortion may not be tolerable. 
The device will not work at low signal levels. 

The receiver blanker is used to eliminate any signals other than 
those received at a preset time for a preset period after the main 
pulse is transmitted. The unit may be switched to use the pulse 
from the tone burst generator, inverted, to cause blanking only during 
the main pulse. It may also be triggered from an externally generated 
delayed pulse in order to discriminate in time against spurious echoes 
in the received signal. Under all circumstances the main pulse must 
be blanked in order to prevent sensitive receiving amplifiers from 
being driven into saturation which can cause unwanted transients for 
some time after the saturating signal is removed. This circuit is 
shown in Fig. 4.6. 

The detector circuit is shown in Fig. 4.7. This unit will 
provide a rapid rise time to follow the leading edge of the echo 
pulse and a long fall time to hold the maximum voltage reached by 
the pulse. The emitter follower input provides the low charging 
resistance with isolation from the input device. The emitter 
follower output provides the high resistance discharge path to hold 
the signal and isolates the detector from the effects of low input 


impedance of any possible following stage. 
91. 
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Fig. 4.9 Calculation of Scattering Area 
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The sampler unit is connected to the detector by a simple SPST 
switch. This allows the detector to be isolated when the sampler is 
used for other purposes such as in the sampling of the wave height 
signal. This unit has been fully described in Section 3 in the 
discussion of the electronics apparatus used for measurement of the 
wave height distribution. The delayed pulse units and the pulse 
height analyzer have also been described in the same section. In 
this application the analyzer measures the distribution of pulse 
voltages and hence gives a measure of the average backscattered 
intensity from the mean and variance of the recorded distribution. 

The regulated power supply for all the supplementary transistor- 
ized gating and sampling circuits is shown in Fig. 4.8. The unit 
incorporates two zener diodes to regulate the voltage at a balanced 


+9v and -9v. The voltage supply is a HP 471A Power Supply. 


Transducers. For the experiments conducted at frequencies 
from 20 to 100 kHz, the projector used was an EDO Model 327 
transducer. This 30 cm diameter transducer contains two concentric 
arrays of ADP crystal piezoelectric elements. When both arrays 
are used simultaneously., the beam width to the 3 db down points at 
180 cm distance is as shown in Table 4.1. The beam angle for these 
3 db down points is also calculated. Similar data are included in 
the table for the transducer used in the experimentat 450 kHz. 
Fig. 4.9 shows the geometry of the calculation of the scattering 
area from the beam measurements. The beam is conical and the area 
of the intersection with the surface is an ellipse. The 3 db 
down points were actually measured for the 70 kHz case and the 


comparison of theoretical and measured data is given in Fig. 4.9. 
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TABLE 4.1 


BEAM ANGLES TO 3 db DOWN POINTS 


Frequency kHz 20 50 70 100 450 
Beam diameter cm 80 43 21 14 10% 
at 180 cm 

Beam angle deg 12.5 6.8 3.3 ታጋ 4,0 


*measured at 72 cm from source. 


In Fig. 4.9 there is a marked departure of the calculated area 
from the measured area. This is due to the effect of spreading 
losses and of the very wide aperture which were neglected in this 
oversimplified model. In the actual calculation of the backscatter- 
ing factor, the measured area was used. Since this area enters the 
expression for the backscattering factor in the denominator, the 
noise line drawn on Figs. 4.17 to 4.21 must be modified to show the 
downward shift the calculation would cause for a constant signal. 

In the experiments at the lower frequencies, it was found that 
the noise level due to interference from the projector power 
amplifier was high enough to mask scattered echoes that could be 
discerned by filtering the broad band signal from a separate hydro- 
phone. Consequently, during the entire run of experiments a separate 
Atlantic Research LC-32 hydrophone, placed at the center of the EDO 
projector, was used. This hydrophone is omnidirectional. The beam 
pattern was modified by placing it against the diaphragm of the EDO 
transducer but in all measurements, this geometry was kept constant. 
Since the backscattering factor is a ratio of scattered to incident 
intensity and both were measured with the same equipment and geometry, 


these effects were effectively cancelled out. 
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For the experiments at 450 kHz, both transducers were 3 cm 
diameter units specially constructed at the Naval Postgraduate 
School. The active elements in each are barium titanate units 
vibrating in a thickness mode. The transducers were mounted side 
by side. As shown in Table 4.1, (at a distance of 72 cm), the 
three db down points on the surface were 10 cm apart making a beam 
angle of 4.0 degrees. Areas calculated from this width and 
modified for the effects of spreading losses at grazing incidence 


were used in the calculation of the backscattering factor. 


Transducer Support and PoSitioning Gear Support and Positioning Gear. The transducers were 


held in specially constructed clamps and secured to a horizontal 
bar which could be accurately (1 degree) positioned to provide any 
angle of incidence for the acoustic beam. The EDO transducer and 
clamp on the bar is shown in Fig. 4.10. The horizontal bar is 
supported by vertical elements which are positioned in height by 
ring clamps. By this means it was possible to maintain the trans- 
ducer at a constant distance from the scattering area. The whole 
was supported by a slotted angle steel bridge which could be 
accurately positioned and leveled at any position over the tank. 


This support bridge is shown in Fig. 4.11. 


Calibration. The need for absolute calibration of the equip- 
ment is avoided by the nature of the backscattering factor defined 
in terms of the ratio of scattered to incident intensity. Only 
the amplifiers in the system with gains to be varied during the 
experiment needed to be calibrated. Particular care was taken to 


ensure that no element in the system was driven into saturation 
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either during the experiment or during the calibration, since 
such an occurrence would cause a change in effective gain even if 


the distortion went unnoticed. 


Comparison With Theory 


Normal Incidence. Equations 9 and 10 in Section 2 show a relation 
between the backscattering factor and the system parameters. The 
roughness parameter (g), is particularly important in these equations. 
The first experiment shows a comparison of this theory which predicts 
the backscattering factor at normal incidence, a specular reflection 
phenomenon, as a function of the roughness parameter. The statistical 
theory was tested for this case by Medwin (15) in 1966. This par- 
ticular experiment compliments that work in that here are statistics 
determined from thousands of echo pulses, and the surface height 
statistics were measured directly. A range of roughness parameter 
from very smooth (g=0.01) to rough (8260) was explored by altering 
the surface roughness with different combinations of blowers and 
by changing the frequency. The theoretical curve and data points 
are shown in Fig. 4.12. The data fit the smooth surface theoretical 
curve closely. The principal source of error is the sensitivity of 
the measurement of the sound reflected from the smooth surface to 
the positioning of the transducer. The actual scattering levels 
were 60 db above the noise. The 450 kHz data at the large rough- 
ness end of the scale are too large. This is considered to be 
caused by the fact that the RMS slope calculated from the optical 
glitter experiment was used to calculate the theoretical curve. Due 


to the large acoustic wavelengths, there are portions of the waves 
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presumably due to capillaries which scatter the light but which 
would not affect the sound. The RMS slope used is too large 


for the acoustic application. 


Backscattered Pulse Distribution. This experiment affords a 
real insight into the nature of the distributions of the scattered 
echoes for specular reflection and for backscatter. Figs. 4.13, 4.14, 
and 4.15 show the echo distributions for normal incidence specular 
scatter at 70 kHz for the rough surfaces with one, two, and three 
blowers in operation. The transition from a normal distribution to 
a Rayleigh distribution as the roughness parameter increases is 
readily apparent. The positions of the 70 kHz data on the curve 
of Fig. 4.12 show that the three different wind conditions give rise 
to roughness parameters occurring over the region of transition from 
smooth scattering to rough scattering. 

Fig. 4.16 is drawn for backscattering data for 20° incidence at 
70 kHz with one blower operating. This is a smooth surface condition 
at 70 kHz, as the specular scatter distribution in Fig. 4.13 shows, 
yet the backscattered distribution is Rayleigh, indicating that the 


backscattered radiation is incoherent. 


Backscattered Pulse Magnitude. Figs. 4.17 to 4.21 show the 


data for backscattering at various angles of incidence under various 
conditions of surface roughness with incident sound projected in 

the upwind direction. On these curves, the predictions of the 
statistical scattering theory are shown. In Figs. 4.17 and 4.18 

the prediction of the Wetzel resonance theory is also shown. This 


curve is not shown in the other figures since the surface wave 
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spectrum was not measured for those surface conditions. The 
Wetzel resonance theory is derived only for relatively smooth 
surfaces. The noise curves are also given to indicate how one 
might inadvertantly measure noise, thinking it to be scattered 
sound, unless particular care is taken to identify the noise level. 
The drop in noise level at high angles of incidence is due to the 
definition of the backscattering factor which is inversely 
proportional to the ensonified area. 

The "tails" of the curves at low grazing angles are observed, 
for the rougher surface conditions, 15 db above the noise. This is 
in an environment free of marine and animal life likely to cause 
volume scattering and free of bubbles by dint of the care taken to 
allow the tanks to settle without agitation or filling, and the 
care taken in filtering to avoid discharge directly into the tank. 
The wind agitation of the surface was not sufficiently violent to 
introduce bubbles into the tank. 

The statistical theory (Section 2, equations 9 and 10) predicts 
the backscattering factor closely at angles near normal incidence 
but fails to predict the leveling of the curves at grazing incidence. 
Although the statistical curves have an inflexion point, it is not 
anywhere near to the actual data above angles of 30 degrees. 

The resonance theory (Section 2, equation 13) using the 
surface wave energy spectrum from Fig. 3.33, has less success in 
predicting the scattering levels in its region of application than 
did the statistical theory in its area of application, yet it is 
clear that the theory resulting from the resonance mechanism of 
backscattering at grazing incidence has the right shape. The theory 


fails the actual measured levels by as much as 10 db, in one case 
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coming too low and in the other case coming too high. The scattering 
predicted by the resonance theory for the predominant frequencies in 
the surface wave spectrum fails to rise to the level of the actual 
scattering in Figs. 4.17 and 4.18 at near normal incidence, indicating 
that the resonance mechanism is unsatisfactory to explain the back- 


scattering in this region. 
Backscattered Continuous Wave Spectra 


The resonance theory of backscattering from a moving surface 
requires the existence of doppler shifts in the scattered radiation 
and further requires that the doppler shifted radiation be due 
principally to certain frequencies in the surface spectrum that would 
give rise to resonances in accordance with Bragg's law. The equations 
for resonance backscattering as developed by Wetzel [9], and modified 
to be consistent with the definition of backscattering factor used 
herein, were used in the last sub-section to predict the backscattering 
near grazing incidence in two of the scattering experiments where the 
data were felt to be sufficiently above the noise level. There it met 
with some success in its predictions. However, it was considered 
that some of the assumptions of the resonance theory should be studied 
experimentally. Some of the results obtained, though largely qualitative, 
yield new insight into scattering analysis and raise some questions 
worthy of consideration in future research in this area. Primarily, 
then, this experiment is directed to three features of scattering: էօ 
find if the doppler shifts can be measured; to find a way to relate the 
spectrum of the continuous wave scattering to the surface wave spectrum; 


to find if there are specific resonances that can be identified. 


109 


ህ:93 5ለ9.- 39299 3ኬ552);1 1३२०२१5 २०४४ Տոօոսլչյսօը 225 "S14 


1524Ճ16ԱՄ 32፲319ሠ7 5973۲1 ۷7۳ 
ƏAPM 1513ԱՅ233340| 


030801103 7٢7 


IL ա - mum | 


| 
| 
| 





— 


1921989813058 
ፈ2399ክከ5932 








55፤3719ሠሃ 
19m0d 





ioonpsuei] 





auoydoapAy 


ከ 
| 
| 
| 


110 


Frequency Synthesizer GR 1161A Coherent Decade Frequency 


Synthesizer 
Power Amplifier HP 467A Power Amplifier 
Transducer EDO Model 327 Transducer 
Attenuator HP 350A Attenuator Set 
Hy drophone Atlantic Research LC-32 
Amplifier HP 463A Precision Amplifier 
Differential Amplifier Tektronix 1A7 in Tek 545 Oscilloscope 
Wave Analyzer GR 1900A Wave Analyzer 


Fig. 4.22(a) Equipment List for 


Continuous Wave Spectra Measurement System 


PT, 


Experimental 


The system used in this experiment is described in the block 
diagram of Fig. 4.22. The frequency synthesizer was stable to better 
than one part in He The wave analyzer local oscillator was 
stable to one part in 10° for short terms of the order of a few 
hours. The CW signal at 48.3 kHz was amplified and supplied to 
the EDO transducer previously described, If there is a doppler shift 
to be detected due to the moving wave system, it should be clearly 
identified by its dependence on the direction of radiation from the 
transducer, whether downwind or upwind. The scattering experiment 
was conducted with the transducer directed for normal incidence 
and subsequently at 60 degrees upwind and downwind to scatter from 
station 3 as shown in Fig. 3.1. The hydrophone was placed approx- 
imately 30 cm from the face of the transducer in clamps to keep the 
distance to the transducer constant. Since the major component of 
the received spectrum would be due to the acoustic radiation incident 
by a direct path on the hydrophone, and since this signal would have 
a constant phase and amplitude relationship to the signal supplied 
to the projector, it could be removed by subtracting a suitably 
attenuated signal derived from the power amplifier, whenever the path 
difference from hydrophone to projector was the correct multiple of 
half-wave lengths to provide the necessary phase relationship. This 
restricted the frequencies that could be used to a set of discrete 
frequencies, but there was sufficient freedom to choose suitable 
frequencies despite this limitation, This incident signal was 
reduced to the point where the signal passing the "skirts" of the 


bandpass filter was much less than any shifted frequency signal in 


MIZ; 


the backscattered spectrum. Since the object of the experiment 
was to look for the doppler shifted frequencies, the loss of 
information concerning any unshifted component of the spectrum 
was of no great consequence. In practice the received signal was 
initially observed under quiet conditions, that is with a smooth 
surface. The difference spectrum under such conditions is shown 
in Fig. 4.23. The specified bandpass characteristic of the wave 
analyzer is superimposed to show that the signal is a single 
frequency so far as can be determined by the analyzer. A second 
component of the signal observed under quiet conditions was the 
reflected signal from the surface. Unless the surface height is 
kept constant during the experiment, the phase relationship of this 
signal to the supplied signal will change and the unshifted 
frequency component measured by the analyzer will vary in amplitude. 
Since the filtered component of the scattered spectrum at any 
frequency is not constant in time but varies widely, the measured 
spectrum is necessarily the result of some sort of averaging process. 
A long time constant filter of the wave analyzer damps the meter 
response. However, the output in this experiment still varied one 
or two db, and the final reading was the result of an averaging by eye. 
In this way, all of the temporal information in the spectrum is 
ignored, yet the process is identical to that which is carried out 
when the surface wave spectrum is analyzed. Hence it is not expected 
that this averaging will invalidate conclusions made about relations 


between the backscattered spectrum and the surface wave spectrum. 
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Fig. 4.24 Difference Spectrum ~ Rough Surface = Normal Incidence 
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Observations 


The result of the spectrum analysis at normal incidence is shown 
in Fig. 4.24. This spectrum is considerably broadened from that 
observed under quiet conditions (Fig. 4.23) and it is generally 
symmetrical, indicating the presence of frequencies shifted from 
the incident acoustic frequency. 

On the other hand, the spectra at 60° incidence both upwind and 
downwind (Figs. 4.25 and 4.26), show decided asymmetry. In the down- 
wind case there is a "bulge" in the negative frequency shifted 
direction; in the upwind direction, the bulge is in the opposite 
direction. The two curves are not mirror images since the down- 
wind spectrum (Fig. 4.26) contains more energy in the shifted 
frequencies than does the upwind spectrum (Fig. 4.25). The slopes 
of these spectra, on the sides not affected by the "bulges", are 
very close to the slope of the spectrum of Fig. 4.24, which was 


obtained at normal incidence, 
Analysis and Speculation 


There is no theory to apply to thís particular phenomenon 
other than general theory concerning surface wave motion, doppler 
shift and frequency domain analysis techniques. In the following 
analysis there are many steps based on conjecture which must be 
examined in greater detail in some future experiment. 

Consider the broadening of the spectrum from the smooth to the 
rough conditions at normal incidence. A number of possible mechanisms 
come to mind. The surface is moving and it should be expected that 


there would be doppler shifts in the echoes from various parts of the 
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surface due to its up and down motion. Since the wave propagation 
is at right angles to the acoustic beam, no doppler shift from this 
motion can be expected. If we consider the received signal in the 
time domain, it is amplitude modulated and, with reference to the 
incoherent nature of the backscatter already considered, it is 
likely that the signal ís also angle modulated. It appears plausible, 
then, that the broadening of the spectrum for normal incidence may 
be interpreted completely as a modulation phenomenon with simultaneous 
amplitude and angle modulation. The spectrum at normal incidence 
was observed to be symmetrical. Since it is possible to generate 
asymmetrical spectra by combined amplitude and angle modulation, but 
not by amplitude modulation alone, it appears that the angle modulation 
must be small. This would be reasonable if the wave height variation 
were small with respect to the acoustic wavelength, a relatively 
smooth surface, but bears examination for rough surfaces. The 
assumption of symmetry of the backscattered spectrum at normal 
incidence is important to the identification of the identification 
of the doppler shifted spectrum. 

Consider the asymmetries which occur in Figs. 4.25 and 4.26. 
It is expected that any doppler shifts due to the moving wave system 
would not be simple since the motion itself is anything but simple. 
However, the wave motion in the rough surface is in one direction 
and the doppler shift must show the correct sign of shift if it is 
to be identified as such. With the wave motion away from the sound 
source, the doppler shift must be negative. In the downwind direction 
the bulge appears on the negative side hence it is an acceptable 


candidate for identification as doppler shift. For sound propagated 
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in the upwind direction, the bulge appears on the positive side; 
further justification for identification as doppler shift. The 
spectra at 60 degrees contain something of the shape of the 
spectrum at normal incidence. This is the broadened spectrum 
about the CW frequency. There is a symmetry that is preserved 
about the CW frequency up to the point where the bulge appears. 
It is postulated, then, that this spectrum is made up of two parts: 
3 symmetrical spectrum due to modulation of the backscattering 
and an additional spectrum due to doppler shift from the moving 
wave structure. It appears to be a reasonable step to attempt 
to separate these two spectra. On the basis of the assumption of 
Symmetry of the backscatter spectrum, the spectral values of the 
unbulged side of the spectrum are subtracted from the values on 
the bulged side. This must be done on a linear scale, not ona 
db scale since the assumption is that the spectra are superimposed. 
Figs. 4.27 and 4.28 show the graphical process of subtraction of 
one sideband from the other (dashed lines) and the resultant 
difference spectra (solid lines). They are ill-defined from 0 to 
6-7 Hz because they are differences of large numbers, and there- 
after appear to have a regular shape not unlike a typical surface 
wave spectrum, However, such an assumption here is totally invalid, 
since the frequency scale is that of the doppler shift. This is 
related to the surface spectrum through the wave velocity, which is 
not related in a simple manner to the surface frequencies. 

In the surface wave velocity curve of Fig. 4.29, it is seen 
that there is a minimum velocity for waves in water due to the 
combined effects of gravity and surface tension on the wave 


formation. Since there is a minimum velocity, there must be a 
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doppler shift frequency. This minimum frequency is readily calculated 
from the minimum velocity of surface waves which is 23.1 cm/sec. 

This corresponds to a surface wave frequency of approximately 12 Hz. 
From Fig. 3.33 it is seen that there is very little energy in the 
surface wave spectrum at this frequency, hence the doppler shifted 
spectrum can be expected to rise from a Zero at or near the calculated 


minimum frequency. The doppler shift equation is 


Af = — (11) 


For a frequency of 48.3 kHz and a minimum velocity of 23.1 cm/sec, 
the minimum doppler shift is 15.1 sin9 where 9 is the angle of 
incidence. 

From the photograph of a typical backscatter at 60" angle օք 
incidence for a 0.2 msec pulse at 48 kHz, (Fig. 4.30), it is seen 
that there are a number of echoes received at times other than the 
time at which the scattered echo in the main beam is expected. The 
time base is 0.5 msec/cm. There is a very large echo starting at 
1.22 ms due to the reflection of the side lobe facing the surface 
at normal incidence. This echo is very wide compared to the 0.2 msec 
pulse width indicating that there are two echoes here. The first at 
normal incidence would produce symmetrical frequency shifts as in 
Fig. 4.24. The other large pulse commences at about 1.3 msec. 
Another pulse is seen at 2.1 msec and the actual pulse from the 
center of the main lobe scattering area occurs at 2.4 to 2.5 msec. 

In the face of this type of scattered return, it is almost impossible 
to make sense of the continuous wave doppler shifted frequencies 
until further consideration is given to the nature of the doppler 


shift. Since the echo at 1.3 msec is very large, the contributions 
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of the pulses at 2.1 and 2.4 msec can be temporarily ignored. 

The effective angle of incidence can then be calculated from the 
times of the return of the normal incidence echo and the major 
scattered echo. Since 1.22 = 1.3 cos9, cos80 = 0.94 from which 

it is seen that the effective angle of incidence is about 20°. 
Using this angle in the doppler shift equation above (1) predicts 
a minimum doppler shift of 5.2 Hz. This compares favorably with 
the minimum frequency of the curves derived in Figs. 4.27 and 
4.28. It must be noted that it is very difficult to determine the 
effective angle of incidence. 

The surface wave spectrum (Fig. 3.33) shows very little energy 
below 1.75 Hz. This indicates that there would be very little 
doppler shifted sound corresponding to velocities associated with 
waves of lower frequency than this. Since the wave velocity of 
gravity waves increases with decrease in frequency, this predicts 
an approximate maximum for the doppler shift spectrum. From 
Fig. 4.29, the velocity corresponding to a frequency of 1.75 Hz 
is 90 cm/sec. This predicts that there will be very little doppler 
shifted sound above 20 Hz. Once again this is in close agreement 
with the maximum observed in the curves constructed in Figs. 4.27 
and 4.28. 

With a certain amount of numerical justification for the 
derived curves, it is possible to speculate on how one might relate 
the doppler shift spectrum, which was conStructed graphically, to 
the surface wave energy spectrum. When equation (11) is combined 
with the equation for the velocity of the surface waves as a 


function of frequency 
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where g is the ecceleration of gravity, v is the velocity of the 
surface wave, and s is the surface tension, a relation is found 
between the doppler shift frequency and the surface wave frequency. 
Due to the implicit nature of (14) this is most easily evaluated 
graphically from Fig. 4.29. By means of this change of abcissa, 

a curve of backscattered doppler shifted energy is plotted against 
the surface wave frequency. This curve is shown in Fig. 4.31 to 
facilitate comparison. The heights of the doppler shifted curves 
for upwind and downwind cases have been set equal to that of the 
surface wave spectrum which is also plotted. The curves fit the 
actual spectrum fairly closely. From the spread of the two 
derived curves it is seen that there is some experimental error. 
The principal source of error is the determination of the effective 
angle of incidence. 

This conjecture seems to go no further, with profit, until the 
assumptions made are subjected to experimental and theoretical 
confirmation. It should be noted that it has been assumed that the 
energy in the doppler shifted spectrum is simply related to the 
energy in the surface spectrum. There is a resonance phenomenon 
which causes the doppler shifted spectrum to alter shape at some 
frequencies. Since there are very few wavelengths of surface waves 


in the scattering area, this resonance has been ignored. 
Conc lusion 


There is a definite asymmetry produced in the backscattered 
spectrum at other than normal incidence, produced by the doppler 


shifts in the backscattered spectrum. An analysis technique has 
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been presented based on assumptions concerning the nature of the 
angle modulation in the backscattered signal and concerning the 
superposition of doppler shifted and symmetrical modulation 

spectra, Ån assumption was implicitly made that the surface 

wave ሟ ቸፍ doppler shifts from each of the spectral components 
moving at its phase velocity. This assumes a linear system for 

the surface waves. The resultant curves of doppler shifted energy 
vs surface wave frequency show some regularity, but there is a 

need for further experimental work and theoretical study before such 
a curve might be used to argue the existence of resonances to 


explain its deviation of shape from the surface waye spectrum. 
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5. Summary 


The phenomenon of backscattering of acoustic radiation from 
a random rough surface has been studied from several points of 
view. Å statistical theory and a resonance theory have been 
examined, and backscattering measurements made to determine whether 
the predicted backscattering factors are close to experimentally 
determined values. The assumptions in theory that the wave heights 
of a rough water surface are Gaussian distributed have been validated 
for this surface as they have been in the ocean itself. In this sense, 
then, the rough surface in the tank is a good model of the sea surface. 
The wave slopes have been shown to be distributed in a Gaussian 
distribution with the upwind-downwind slope being skewed in the down- 
wind direction. The upwind-downwind to crosswind slope ratios are 
slightly higher than those experienced at sea when the rougher surfaces 
are used. This is a measure of the effect of creating the rough 
surface in a channel. 

Over a very wide range of roughnesses (0.01 <g<60.0), the 
specular scattering at normal incidence is closely predicted by a 
statistical theory due to Beckmann and Spizzichino (5). The 
statistical theory predicts near-normal backscatter reasonably 
well, but it fails completely to predict backscatter for incidence 
angles greater than 30°. The predictions concerning the distributions 
of scattered pulse heights have been substantiated by actual measure- 
ments: the near coherent scattering is Gaussian distributed and the 
incoherent is Rayleigh distributed. 

At grazing incidence, the shape of the backscattering factor 


curve is fitted closely by the resonance theory prediction. There is 
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some discrepancy ín the predícted levels which may be due to 
inaccuracies in measuring surface wave spectrum and in the measure- 
ment of backscattering factor from experimental data. A theory of 
backscattering would do well to consider both mechanisms for 
scattering as have been studied. Each mechanism has a dominant 
effect in its own region of application, and the actual backscatter- 
ing phenomenon is closely predicted by a summation of backscattering 
factors from the two theories. 

Further consideration of backscattered signals as a modulation 
phenomenon may yield a method of extracting the surface wave spectrum 
from the backscattered spectrum. An analytical method has been 
suggested, but it is based on several assumptions which have not yet 


been validated from theory or in experiment. 
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